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2 ■ ABSTRACT 
> . 

\ Azimuthal color (age) gradients across spiral arms are one of the main predic- 

^ \ tions of density wave theory; gradients are the result of star formation triggering 

\ by the spiral waves. In a sample of 13 spiral galaxies of types A and AB, we 

^ \ find that 10 of them present regions that match the theoretical predictions. By 

\ comparing the observed gradients with stellar population synthesis models, the 

I pattern speed and the location of major resonances have been determined. The 

00 \ resonance positions inferred from this analysis indicate that 9 of the objects have 

^ \ spiral arms that extend to the outer Lindblad resonance (OLR); for one of the 

\ galaxies, the spiral arms reach the corotation radius. The effects of dust, and 

^ \ of stellar densities, velocities, and metallicities on the color gradients are also 

c3 ' discussed. 
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Introduction. 



Density A vave piienomena iiave been proposed to expla i n tiie spiral structu re seen in 
disk galaxies JLindbladlllQGsl : llin fc Sliulll964l : [Toomrelll977l : ISertin et aLlll989al lbll Obser- 
vationally, these phenomena are best studied in the nea r infrared (near-IR ), especially the 



i^'-band, that mostly traces the old stars in the disk (e.g.. lRix &: Riekdll993l ). The ol d stellar 
disk, however, has a disordered optical counterpart of young stars, gas, and dust (IZwicky 
19551 : iBlock fc Wainscoatlll99ll ). so the question arises whether the spiral structure, seen in 
near-IR bands, and the star formation, seen in the optical bands, are coupled or not. If disk 
dynamics and star formation are indeed related, the star formation rate per unit gas mass 
should be affected by the presence of the density wave. Unfortunately, H2 cannot be directly 
observed, and the r elation betw een detected CO and total H2 mass is a whole controversial 
topic in itself (e.g 



Allenlll996l . and references therein). 



The alternative approach of comparing star formation rates, past (as traced by optical 
and near-IR surface photometry) and present (as probed by Ha emission), in galaxies with 
different Hubble types has been carried out; as a result, both a positive correlation between 
disk dynamics a nd star formation (e.g . . ISeigar fc Jamej 120021). and t he absence of such a 
correlation (e.g.. lRyder fc Dopitalll994l : lElmegreen fc Elmegreerull986l ) have been claimed. 



In this work, we will focus on the relation between star forrnation and density wav e phe- 
nomena, as proposed in the large scale shock scenario ( Roberts 1969 : Shu et aP 1972). Evi- 



dence has been ga thered, from observation s of dust, gas compression (e.g.. 



Mathewson et al. 



I972I : IVisserl Il980l ) , and molecular clouds (IVogell Il988l : ISchinnerer et al.l 120041 ) near the con- 
cave regions of spiral arms, that suggests that star formation is triggered there. Through a 
simple interpretation of these ideas and observations, the existence of azimuthal color-age 
gradients has been predicted. As seen in Fig. [H if we assume that the spiral pattern rotates 
with constant angular speed, and that the gas and stars have differential rotation, a coro- 
tation region exists where these two angular speeds are equal. At smaller radii, bursts of 
star formation take place where the di fferentially rotating gas overtakes the spiral wave, an d 
show up as brillant HII regions (e.g., [Morgan et al.lll952l : lElmegreen fc Elmegreenlll983al ). 
As young stars age, they drift away from their birth site, thus creating a color gradient. Star 
formation can occur also beyond the corotation radius, when the spiral pattern overtakes 
the gas. The assumption that the patt ern rotates with constant angular speed can be cor- 



roborated w ith numerical simulations (IThomasson et al.l Il990l : iDonner fc ThomassonI Il994 



Zhang||l998l ). but the only observational evidence of this premise is the apparent persistence 



of the spiral structure for up to a Hubble time in nearby galaxies ( lElmegreen fc Elmegreen 
1983bl ). in avoidance of the winding dilemma. 



Another way to test the assumption of constant angular speed of the spiral pattern 
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is by the dynamical consequences it has on the disk environment and material. In this 
regard, the prediction of azimuthal color gradients can be seen as another test of the spiral 
density wave theor y. So far, the s earch for these color gradients has in general yielded 



1990 



inconclusive results (ISchweizerlll976l : iTalbot et al.lll979l : ICepa fc Beckmanlll990l:lHodge et al. 



del Rio fc Cepa 19981 ). and a few e xceptional positive ca ses hke lEfremovl msS ) for 



M31, ISitnikI (119891 ) for the Milky Way, and lGonzalez fc GrahamI (119961 ) for the spiral galaxy 
M 99. 



Gonzalez &: GrahamI (1l996l . GG96 hereafter) used for the first time a supergiant sensi- 



tive and reddening-free photometric index0 to trace star formation, defined as 



Q{rJgz) = {r-J)-^^^{g-z) 



Using the extinction curves of ISchneider et al.l (119831 ). and iRieke fc LebofskyI (119851 ) for a 
foreground screen, the color excess term is E{r — J)/E{g — i) = 0.82. According to pop- 
ulation synthesis models, following a star formation burst the Q{rJgi) photometric index 
increases its value for ~ 2.6 x 10^ years, and then starts to decline, as shown in Fig. [2] for 
most recent Chariot & Bruzual (2007, in preparation, CB07 hereafter) models with different 
stellar mixtures in which young stars constitute from 1% to 5% by mass. Dust lanes arising 
from spiral shocks are expected to precede the azimuthal color gradients, hence the impor- 
tance of reddening- f ree di agnostics. With the stellar population synthesis (^SPS) models of 



Bruzual fc CharlotI (I19931). a nd the radiative transfer models of iBruzual. Magris fc Calvet 



mm and IWitt et all (jl992|), GG96 investigated the behavior of Q{RJVI) for mixtures of 
dust and stars with different relative spatial distributions. The main effect of a mixture 
of dust and stars on the Q index, illustrated in Fig. [3] for the CB07 population synthesis 
models, is to increase its value relative to the r = and foreground dust screen cases, and 
Q{RJVI) is no longer reddening-free when ry > 2§| H owever, the requirement that Ty < 2 



is lik e ly fulfilled by the di sks of nearly face-on galaxies (jPeletier et al.lll995l : iKuchinski et al. 



1998 



[e.g., 



Xilouris et al 



Chariot fc Fall 



9991) : moreover, as stars age, dust dissipates, and reddening diminishes 



2000) 



Based on their study of M 99, GG96 propose that real data are best matched by SPS 
models with 0.5% to 2%, by mass, of young stars. These values are in accordance with 



^ See Binney & Merrifield (1998), for other reddening- free indices. 

^The same conclusion can be applied to QirJgi), since the optical bands R, V, and / have approximately 
the same effective wavelengths, respectively, as the r, g, and i filters; the main difference between both sets 
is that the latter passbands are narrower. 
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Fig. 1. — Stellar age gradients across the spiral arms are indicated by arrows that go from 
blue to red. The azimuthal age gradients are produced by stars born in the spiral shock, 
where the shocked interstellar medium forms a dust lane, that later drift away as they age. 
The direction of the gradients changes at the corotation radius, -Rcr- Inside this radius, the 
disk material overtakes the spiral wave, and beyond i t the spiral wave catches up with the 
material [see also figure 1 in lPuerari fc Dottoril (119971 )]. 



SxlO' 10» 1.5x10" 

Fig. 2. — Theoretical Q{rJgi) vs. time, CB07 models. The duration of the burst is 2 x 10^ 
years, with a Salpeter IMF; the fraction of young stars ranges from 1% to 5% by mass. Lower 
and upper mass limits are 0.1 and 10 Mq, respectively. 
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Fig. 3. — Model Johnson fil ter Q(RJVI) vs. t ime (t) for CB07 models, reddened as per the 
"starburst galaxy" model of IWitt et al.l (1l992l ). The duration of the burst is 2 x 10^ years, 
with a Salpeter IMF and 2% by mass of young stars. Lower and upper mass limits are 0.1 
and 100 Mq, respectively. 
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the am ount of B light in the arms contributed by young stars, as estimated by ISchweizer 
(Il976l ). GG96 also hypothesize there is an inverse correlation between high-mass star- forming 
regions and detectable azimuthal color gradients. The color gradient in M 99 lies where no 
HII regions are identified in Ha images of the galaxy. This counterintuitive finding could help 
explain the dearth of positive detections to date, if contamination from bright emission lines 
produced in HII regions around the most mass ive stars masks the color gradients associated 
with star formation in spiral shocks (jShulll997l ). On the other hand, evidence has been found 
recently of star forming regions w ith very few massive stars that generate only scant Ha 
emission ( llndebetouw et al.ll2008l ). In this investigation we apply the GG96 method, using 
the Q{rJgi) photometric index, to a new sample of galaxies in order to search for and analize 
color-age gradients near spiral arms. 



1.1. Color gradients: the link between star formation and spiral dynamics. 



(Lin 


1970; 


Mark 


1976; 


Lin & Lau 


1979; 


Toomre 


1981; 



Of these resonances, the most important ones are the inner Lindblad resonance (ILR), the 
4:1 resonance, corotation (CR), and the outer Lindblad resonance (OLR)i 

Azimuthal color-age gradients retain some information about the stellar drift relative to the 
spiral shock, allowing us to obtain Qp, the angular velocity of the spiral pattern. In order to 
find Qp from the gradient information, we can use: 



M^cot(-,)ln(||). (3) 

where t is the age of the young stellar population at an angle 6'shock away from the shock 
position (t' is a variable of integration), v{t') is the velocity vector of the young stellar 
population in an inertial reference system; (p{t') is the unit vector in plane polar coordinates 
p, (f, in a non inertial reference system; R{t') is the orbital radius of the studied regioii, 
measured from the center of the disk to the center of mass of the young stellar populationp 

At the ILR (OLR), the epicychc frequency k = ±2(il — k = 4(f7 — ftp) at the 4:1 resonance. 
'^Rit) corresponds to 6'shock- 
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and i is the arm pitch ang lei if the azimuthal angle increases in the direction of rotation 
and the spiral arms trail. The angular quantity accounts for the logarithmic spiral shape 
of the shock. Assuming the departures from circular motion are small, we have: 

U^.^Pl^A, (4) 

where Rmean is the mean orbital radius of the studied region, Vrot is the circular orbital 
velocity in an inertial reference system, and t is the age of the young stellar population at 
the azimuthal distance d from the shock. Expanding the term Vrotit')dt', in a Taylor series, 
we get : 

O _ 1 f iVrot)\ot + U'-^)\ot' + ( 
m.p.n.n. \ ^ 



where C represents hi gher order terms in the expansion. According to theory (e.g., [Roberts 



19691 : ISlyz et al.ll2003l ). the higher order terms may account for 20-30 km s~^, a quantity that 
is of the same order as the mean rotation velocity error, after considering the uncertainty 
due to galaxy inclination (see Table [3]). Hence, we can neglect the higher order terms in 
eq. [HI such that Vrot ~ constant, and obtain 



d 



= ( Vrot --]. (6) 



VLp is found by stretching the model (which gives Q as a function of t) to fit the data (where 
Q is a function of d). 



2. Observations and data reduction. 



Our total sample of objects consists of 31 almost face-on spiral galaxies of various 
Hubble types, with angular diameters between A' and 6^ T his sample was chosen from 
the Uppsala general catalogue of g alaxies (UGC; [Nilson|[l973l ). the ESO-Uppsala survey of 
the ESO (B) atlas (lLaubertsl[l982l ). and the Second reference catalogue of bright galaxies 



''The angle between a tangent to the spiral arm at a certain point and a circle, whose center coincides 
with the galaxy's, crossing the same point. 
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( jde Vaucouleurs et al.l Il976l ) . From this sample, w e further select for their analysis in the 



present paper 13 A and AB (jde Vaucouleurs 



19591 ) galaxies, based on the visual inspection 



of their 2-D images in the Q{rJgi) photometric index. In this diagnostic, the disks of some 
of the galaxies in our original sample of 31 clearly appear divided in two halves, each one 
with a different average value of Q. The analysis of these "Q effect" galaxies and of the 
remaining barred galaxies will be undertaken in subsequent publications. 

The data were acquired during 1992-1995 with four different telescopes: the Lick Ob- 
servatory 1-m, the Kitt Peak National Observatory (KPNO) 1.3-m, and the Cerro Tololo 
Interamerican Observatory (CTIO) 0.9-m and 1.5-m telescopes. De ep photometric ima ges 
wer e taken in the optical filter s g, r, and i, and in the near-IR J, Kg (jPersson et al.l 1 19981 ) or 
K' (IWainscoat &: Cowielll992l ). Effective wavelengths and widths of all the filters are listed 
in Table (H the observation log for the 13 galaxies is shown in Table [21 

The CCD at the Lick 1-m telescope was a Ford 2048^, with a pixel scale of (y.'185 pixel"^. 
For the infrared observations at Lick, the same telescope was fitted with the LIRC-2 camera; 
it had a 256^ NICMOS II detector, with a 1'.'145 pixel"^ plate scale. The CTIO 0.9-m 
optical telescope used a Tek 1024^ and a Tek 2048^ CCDs, both with a 0'.'4 pixel"^ plate 
scale. The CTIO infrared observations were performed at the 1.5-m telescope, with the 
CIRIM instrument, which used a 256^ NICM0S3 array; the CIRIM focus was adjusted to 
give a I'.'IG pixel"^ plate scale. The KPNO infrared observations were made with the IRIM 
camera, that employed a 256^ NICM0S3 array, with a 2" pixel"^ plate scale. 

The data were reduced with the image processing package IRAI^ ( Tody 1986 . 1993), 
using standard techniques. For the optical data, overscan and trimming corrections were 
first applied. Bias subtraction and fiat field division corrections were used. Master fiats were 



^ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 



Table 1. Filter characteristics 



Filter 




FWHM 


9 


5000A 


830A 


r 


6800A 


1330A 


i 


7800A 


1420A 


J 


1.25/im 


0.29/im 


Ks 


2.16/im 


0.33/im 


K' 


2.11ftm 


0.35/^m 
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produced from stacks of twilight flats that were averaged pixel by pixel, after scaling each flat 
by its median value and sigma-clipping deviant pixels; the master flats were then divided by 
their mean pixel value. For sky subtraction, the sky level was determined by masking bright 
objects in the image and subsequently fitting a constant value to the remaining pixels, after 
rejection iterations. To produce mosaics, individual frames were superpixelated by a factor 
of 2 in each dimension, and registered to the nearest half (original) pixel. Before adding 
into the final mosaic, cosmic rays were removed by comparing each individual image with a 
median one. 

A correction for non-linearity was applied to the near-IR data. The correction was 
obtained by adjusting a polynomial function, pixel by pixel, to dome fiats of increasing 
exposure times, taking into account variations in the count-rate between exposures. During 
observations, frames were mostly taken in the sequence SKY-SKY-OBJECT-OBJECT-SKY- 
SKY. . . Objects in the sky were masked, and the sky frames were median scaled, averaged, 
and subtracted from the object. This procedure also takes care of dark current removal. Flat 
field corrections were applied with master fiats derived from dome fiats. Cosmic ray removal 
and mosaic registering were done with the same procedure used for the optical data. 

Th e optical calibration was performed usiri g synthetic photometry in the Thuan-Gunn 
system ( iThuan fc Gunnlll976l : IWade et al.lll979l ). This system is based on the standard star 
BD-|-17°4708, to which the magnitude g = 9.5, and the color indices g — r = r — i = are 
assigned. Synthetic mag nitude; for other spectrophotometric ^t^dardfl were obtained by 
means of the following relation: 



mag, 



syn 



-2.5 log 



10 



(7) 



where f\ is the absolute spectral energy distribution of the star in ergs s ^ cm" -2 A" ; ZP IS 
the zero point, and -R(A) is the system response curve, including the quantum efficiency of the 
detector, the transmission function of the filter, and the effect of atmospheric absorption. 
The atmospheric t ransmission curve for the northern h emisp here observations was taken 
from Hayes ( 1970 ). and for C TIO from iHamuv et al" (1992). The absolute fluxes were 
obtai ned, for BD+17°4708 fromlOke fc GunnI (|l983r). and for t he re m aining standa r d star s 
fro mlstonel Jl977 ) iMassev et all (|l988l ). iMassev fc Gronwalll (|l99oh . Iflamuv et all (|l992l ). 
and lHamuy et al.l (Il994j ). Individual galaxy frames taken under photometric conditions were 
used for calibration, and mosaics were scaled to such photometric data. 



^ Feige 15, 25, 34, 56, 92, 98; Kopff 27; LTT 377, 7987, 9239; EG 21; BD+40°4032; and Hiltner 600. 
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For the infrared data, we only had one photometric season at KPNO and none at 



CTIO. In order to obtain a uniform photometry, we ca 
images from the Two Micron All Sky Survey (2 MASS, 



ibrated our J and data with 



Skrutskie et al.l Il997l. mm. For 



the K' data, we use photometric standard s from iHawarden et al.l (I2OOII ). and adopt K' = 
K + 0.2{H — K) (jWainscoat fc Cowielll992l ). We do not include color correction terms in our 
infrared calibration, but take this systematic error into account in the zero point uncertainty. 



Finally, the optical images were degraded to the lower resolution of the infrared images, 
and aligned with them in order to proceed. 
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Table 2. Observation Log 



Object 



Filter Exposure(s) 



Telescope 



Date (month/year) 



NGC 4939 



NGC 3938 



NGC 4254 



NGC 7126 



NGC 1417 , 



NGC 7753 



NGC 6951 



NGC 5371 



9 


2700. 


CTIO 0.9 m 


3/94 


g 


300. 


Lick 1 m 


4/94 


r 


2400. 


CTIO 0.9 m 


3/94 


r 


300. 


Licit 1 m 


4/94 


i 


3600. 


CTIO 0.9 m 


3/94 


i 


300. 


Lick 1 m 


4/94 


J 


840. 


Kitt Peak 1.3 m 


3/94 


Ks 


515. 




55 


9 


600. 


Lick 1 m 


4/94 


r 


300. 




" 


i 


1189. 




2/94, 4/94 


J 


1252. 


Lick 1 m 


2/95 


J 


600. 


Kitt Peak 1.3 m 


3/94 


Ks 


560. 


" 


" 


g 


4443. 


Lick 1 m 


4/94 


r 


2100. 


5) 


55 


i 


3308. 




55 


J 


501. 


Lick 1 m 


2/95 


J 


960. 


Kitt Peak 1.3 m 


3/94, 11/94 


Ks 


590. 




55 


9 


3900. 


CTIO 0.9 m 


9/94 


r 


3900. 


5) 


55 


i 


3900. 




" 


J 


300. 


CTIO 1.5 m 


9/95 


Ks 


390. 


n 


9/95 


9 


3600. 


Lick 1 m 


9/93, 10/93 


r 


5048. 




9/93, 10/93, 11/93 


i 


3600. 


" 


9/93, 10/93 


J 


720. 


Kitt Peak 1.3 m 


9/93 


Ks 


256. 






9 


4500. 


Lick 1 m 


11/93, 10/94, 11/94 


r 


5812. 


)) 


55 


i 


5307. 




55 


J 


626. 


Lick 1 m 


10/94 


K' 


113.4 






9 


6120. 


Lick 1 m 


6/92, 8/92, 9/93 


r 


6900. 






i 


6296. 


5) 


7/92, 8/92, 9/93 


J 


2136. 


Lick 1 m 


9/93, 11/94 


J 


4140. 


Kitt Peak 1.3 m 


7/94 


Ks 


608. 


5) 


9/93, 11/94 


9 


5100. 


Lick 1 m 


6/92, 3/93, 4/94 


r 


6600. 


)! 


55 


i 


2375. 


5) 


6/92, 4/94 


J 


501. 


Lick 1 m 


2/95 


J 


870. 


Kitt Peak 1.3 m 


3/94 
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Prior to analysis proper, the images were deprojected using position angles (P-A.) and 
isophotal diameter ratios ta ken from the Third Reference Catalogue of Bright Galaxies (RC3, 
de Vaucouleurs et al.lll99ll ). unless indicated otherwise in Table O The inclination angle (a) 
was obtained using the simple approximation cos a = [b/a). 



The spiral arms where t hen unwrapped by plotting them in a ^ vs. In R map (jlye et al. 



19821 : lElmegreen et al.lll992n . Under this geometric transformation, logarithmic spirals ap- 
pear as straight lines with slope = cot{-i), where i is the arm pitch angle. Following the 
procedure used by GG96, the phase of 9 was then changed for each (fixed) value of In R, until 
the arms appeared horizontal. This way, selected regions with candidate color gradients can 
be easily collapsed in In R to yield 1-D plots of intensity vs. distance; features in such plots 
have a higher signal-to-noise ratio than in 2-D images, and can be directly compared to the 
SPS models. 



3. The Q{rJgi) photometric index and the stellar population synthesis models. 



With the aim of tracing star formation across the spiral arms of disk galaxies we use 
the Q{rJgi) reddening- free photometric index, as defined by GG96. The Q{rJgi) index was 
chosen between different possible filter combinations for two reasons: (1) the quotient of 
relevant color excesses is close to un ity, and hence all bands involved have similar weights, 
and (2) because the SPS models (e.g.. lBruzual fc Charlotll2003l ) predict a detectable gradient. 
To explain how Q works, GG96 express it in logarithmic form: 



j2.05j2.50 

QirJgt) = logio p.sorlos ' (8) 

r i 

where Ix is the light intensity in each filter. In active star forming regions, blue and red 
supergiants dominate the scene, and the Q{rJgi) index has higher values because the g 
and J bands in the numerator of eq. [H] are tracing, respectively, the light from these stars. 
Conversely, when star formation activity is poor, the Q{rJgi) index has a relative lower 
value. 

For the present work, we use a preliminary version of the CB07 models for comparison 
with the observations!^ A star formation burst was added to a background population of old 



^ The C jB07 code uses the Padova 1994 single stellar population evolutionary tracks as assembled and de- 
scribed by iBruzual fc Charlotjboosh . but includes the new prescription for s tellar evolution in the thermally- 
pulsating asymptotic giant branch (TP-AGB) bv iMarigo fc Girardil (|2007r ). Some differences between the 
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Table 2 — Continued 



Object Filter Exposure(s) Telescope Date (month/yeax) 



Ks 520. 

NGC 3162 g 5400. Lick 1 m 3/93, 11/93, 4/94, 10/94, 11/94 

r 6000. " 4/93, 11/93, 4/94, 10/94, 11/94 

i 1800. " 11/93, 4/94, 10/94, 11/94 

J 990. Kitt Peak 1.3 m 3/94 

Ks 560. 

NGC 1421 g 3600. CTIO 0.9 m 9/94 

r 3600. 

i 3550. 

J 2241. Lick 1 m 10/94, 12/94 

J 330. CTIO 1.5 m 9/94 

J 300. Kitt Peak 1.3 m 11/94 

Ks 160. CTIO 1.5 m 9/94 

Ks 60. Kitt Peak 1.3 m 11/94 

NGC 7125 g 3900. CTIO 0.9 m 9/94 

r 3900. 

i 3900. 

J 600. CTIO 1.5 m 9/95 

Ks 285. 

NGC 918 g 4500. Lick 1 m 11/93, 10/94, 11/94 

r 4500. 

i 5400. 

J 2567. Lick 1 m 10/94, 12/94 

K' 897.6 

NGC 578 g 3900. CTIO 0.9 m 9/94 

r 3600. 

i 4800. 

J 240. CTIO 1.5 m 9/95 

Ks 375. 
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Table 3. Galaxy parameters 



Name 


Type 


P.A. (degrees) 


ajh 


Vma-x (km S ^) 


Radial velocity (km s ^) 


Distance (Mpc) 


NGC 4939 


SA(s)bc 


10 


1.95 ± 0.13 


207.5 ±7.6 


3111 ± 5 


46.5 ±4.0 


NGC 3938 


SA(s)c 


52=' 


1.10 ±0.05 


39.2 ±3.6 


809 ±4 


15.8 ±1.4 


NGC 4254 


SA(s)c 


68*= 


1.15 ±0.05 


106.7 ±6.4 


2407 ± 3 


16.5 ±1.1== 


NGC 7126 


SA(rs)c 


80 


2.19 ±0.15 


152.4 ±2.8 


3054 ± 8 


44.6 ± 3.8 


NGC 1417 


SAB{rs)b 


10 


1.62 ±0.11 


205.1 ±4.3 


4120 ± 17 


57.1 ±4.8 


NGC 7753 


SAB(rs)bc 


50 


1.58 ± 0.15 


164.8 ±9.5 


5163 ± 4 


72.1 ±6.1 


NGC 6951 


SAB(rs)bc 


170 


1.20 ±0.08 


144.2 ±6.0 


1426 ± 6 


24.9 ±2.1 


NGC 5371 


SAB(rs)bc 


8 


1.26 ±0.09 


177.0 ± 11.8 


2553 ± 7 


43.5 ±3.7 


NGC 3162 


SAB(rs)bc 


31'' 


1.20 ±0.08 


76.6 ±7.8 


1298 ± 7 


23.7±2.1 


NGC 1421 


SAB(rs)bc 





4.07 ±0.28 


165.9 ±8.8 


2090 ± 5 


29.3 ±2.5 


NGC 7125 


SAB{rs)c 


110 


1.45 ± 0.10 


104.2 ± 18.2 


3054 ± 8 


44.6 ± 3.8 


NGC 918 


SAB{rs)c 


158 


1.70 ±0.16 


112.5 ±6.2 


1509 ± 4 


21.7 ± 1.8 


NGC 578 


SAB{rs)c 


110 


1.58 ± 0.07 


117.8 ±4.1 


1630 ± 4 


22.7 ± 1.9 



Note. — Col. (2) and (3). Types and position angles from RC3. Col. (4). Isophot al diameter ra t io der ived from the iJ25 
parameter in RC3. Col. (5). Maximum rotation velocity obtained from the HI data of lPaturel et al.l ||2003|) . not corrected for 
inclination. Col. (6). Heliocen tric radial velocity from RC3. Col. (7). Hubble distance obtained from the heliocentric radial 
velocity and the infall model of lMould et al.l l l2000h . 

HPaturel et"al] ll200d') 

IPhookun et al.l l ll993h 

'^Distance to Virgo from lMei et al. I liool), adopted since NGC 4254 (M 99) is a member of the Virgo cluster (the procedure 
used for the other galaxies yields a Hubble distance of 40.2 ± 3.5 Mpc). 
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stars with an age of 5 x 10^ years; both components have a Salpeter IMF. Using an older 
background population makes no difference in Q{rJgi). The burst has a duration of 2 x 10^ 
years, in accordance with the age spread of OB associations inferred from observations (e -g-i 
Elme green fc Ladal 119771 : IPoom et al.lll985l : iMassey et al.lll989l : iBurningham et al.l 120051 )1^1 



Each of the four bands needed to produce the model Q{rJgi) index was calculated as 
shown below for g: 



g,,p{t) = -2.5 logio(/3/10-°-^^--W + /3,,10-°-^^°'^), (9) 

where gcspit) is the ^f-band surface brightness, in mag arcsec"^, of the composite model 
stellar population (csp) made up of young and old stars, at time t after the formation of 
the young stars. Pi and Pu are, respectively, the fractions of young and old stars by mass 
{Pi I = 1 — Pi); gyngit) is the g surface brightness of the young single stellar population at 
time t] and (^oid is the (/-band surface brightness of the old single stellar population. The 
same holds for the r, i, and J bands. 



4. Results. 



For the time being, we adopt models with constant stellar velocity (i.e., stellar orbits 
are assumed to be circular) and with a constant fraction of young stars of 2% by mass. IMF 
upper mass limits of both 10 and 100 Mq are considered. No data were matched by the 
models with Mupper = IOOMq, however. We also use models with solar metallicity in both 
populations, and shift them to the surface brightness level observed from the data for each 
object. A possible reason for the need of this shift is metallicity, as will be explained in § 15.31 



two sets of models are illustrated in lEminian et al. ( 20081 ). 



^The star formation timescale in these OB associations is a raa-tter 
netic field-regulated model of star formation (jShu. Adams fc Lizano 

t he range 5-10 Myr, while in the supersonic turbulence-regulated star formation model (e.g., 



of debate. The mag- 

19871) predicts timescales in 



Ballesteros-Paredes. Hartmann fc Vazauez-Seinadenilll99 9') the process up to the of pre-main-sequence stars 
takes around 3 Myr. Mouschovias et al. ( 2006h discuss observational evidence from external galaxies, indi- 



cating that the lifetime of molecular clo uds and the timescale of star formation are ~ 10'' years. However, 



Ballesteros-Paredes fc Hartmanij (|2007l ) propose that the large age spread observed in OB associations is 
due to successive generations of stars, triggered by the effects of stellar energy input through photoion- 
ization, stellar winds, and supernovae. This process eventually slows down as gas gets dispersed and the 
gravitational potential changes. Our adopted time for star formation (2 x lO*" years) is in accordance with 
the observational results, although the details of the star formation processes involved may differ. 
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The rigorous application of more complex models would require an individual treatment of 
each galaxy and region. 

Once Qp is known, it is possible to locate major resonances in the galaxy deprojected 
images, thus providing a link between the color gradients and the dynamics of the spiral 
disk. That is, if the gradients are indeed caused by star formation in large-scale shocks, the 
positions of the spiral endpoints must be consistent with the location of major resonances, 
as deduced from a comparison between the derived flp and the known orbital velocity. 

There are rotation curves available in the literature for a few objects in our sample, but 
observed w ith different cr i teria. With the purpose of homogenizing our sample, we use the 



HI data of iPaturel et al.l (120031 ). who give the maximum rotation velocity uncorrected for 
inclination. These values are shown in Table [31 An important aspect to notice in equation [H] 
is that Rmean, Vrot, and d depend on the inclination angle a used to deproject the images. 
An independent variable form of this equation is shown in the Appendix (see equation lAip . 
which we use to obtain the uncertainty in Qp. Also, in order to d e termi ne Qp, the distance 



to the object must be known. We use the model of iMould et al.l (120001 ) to convert galaxy 



heliocentric velocities, obtained from the RC3 catalogue, to Hub ble flow velocities; from 



1 

these we find the distance, adopting Hq = 71 ± 6 km s~^ Mpc~^ (IMould et al.ll2000l ). At 
small redshift, the uncertainty in the model is expected to be ~ 100 km s~^. The uncertainty 
for any single galaxy is about 100 km s~\ plus the field dispersion (perhaps as big as 250 
km s~^, and much larger in cluster cores), added in quadrature. The model works well for 
groups and clusters locally, and less well (by definition) for individual galaxies, unless they 
are in a quiet part of the flow (Huchra, J. P. 2008, private communication). The heliocentric 
velocities for our objects and the distances derived from this model are shown in Table [31 
To obtain the distance uncertai nty we propa gate the errors of the heliocentric and infall 



velocities in the equations of the [Mould et al.l model. 



All the spiral arms visible in the mosaics were inspected for profiles similar to those 
expected from azimuthal color gradients. Only a few regions across the arms of the analyzed 
objects present profiles that match the theoretical expectations. These regions are marked in 
the optical mosaics of each individual object. The Q{rJgi) index data, tracing star formation; 
the {g — J) data, outlining the dust lane location, and the Kg (or K') data, following the 
density wave, are shown for each single region. The site with the highest {g — J) value (i.e., 
the highest extinction) is taken as the origin of the (azimuthal) distance (in kpc), which 
increases in the sense of rotation; the direction of rotation was inferred assuming that arms 
are trailing. The comparison between model and data Q{rJgi) index is shown in separate 
figures, where the vertical error bars correspond to ±lcr of the computed uncertainty in Q, 
including read-out noise and sky subtraction. The horizontal bars represent the uncertainty 
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in the "stretch" apphed to the the model in order to fit the data. This is ultimately an 
uncertainty in the stellar drift velocity and, hence, in Qp, and is obtained as explained in 
the Appendix. We include in this error an estimate of the different stretch that would be 
required by models with variable densities, velocities and metallicities discussed in § [S] (see 
figure [TTj) . The resonance locations determined once Qp was obtained from stretching the 
stellar population model to the data are marked on the infrared mosaics. The observed 
parameters for each region and the derived dynamic ones are summarized in Table HI The 
Q index magnitude offset, applied to the solar metallicity model with 2% of young stars in 
order to fit the data, is also listed. 

From our Kg and K' data we obtain visually the location of the spiral endpoints in our 
objects. In some few cases the signal-to-noise ratio was too low, and the J or the optical 
i images were used instead. The positions of the derived spiral endpoints, as well as the 
wavelength used, are listed in Table |H We compare these spiral endpoints with the locations 
of major resonances (4:1 resonance, corotation, and OLR) inferred from the Q{rJgi) index 
data and the stellar models, as described in § 11.11 For objects with more than one studied 
region we choose the one with the lowest uncertainty and the best match with the spiral 
endpoints. The selected regions are marked with an asterisk in Table |H The results of this 
analysis are summarized in Fig.|H the vertical axis corresponds to the Rres/Rend ratio, where 
Rres is the location of the major resonance derived from the analysis, and Rend is the spiral 
endpoint obtained from the data. 

With the exception of NGC4254 (M 99), objects are organized by Hubble type (from 
SAbc to SABc). Remarks in the caption of Fig. [18] apply also to figures [T9] - |4T| unless 
indicated in individual captions. 

NGC 4254 : [Figures M - [IS] Although NGC 4254 is a "Q effect" galaxy with two 
asymmetric halves in this diagnostic, we include it here as a test of the consistency of 
our procedure vis-a-vis its first application by GG96. The resonance positions computed 
for region NGC 4254 A (the same analyzed by GG96) place the spiral endpoint for the 
corresponding arm at the OLR. Even though the Qp for region NGC 4254 B differs from 
the value obtained for region A, the spiral endpoint of the arm to which this region belongs 
matches the position of its OLR tooO Not surprisingly for a "Q effect" galaxy, region B 
Q values are also 0.03 magnitudes below those for region A. Given the fact that the two 
gradients yield different values for Qp but their corresponding arms end at their respective 



^°Results differ wfien adopting Vrot = 140 km s ^ and a — 42°, as derived by Ptiookun et al. (1993) and 
used by GG96. Witli ttiese values, ftp = 19.0 ± 2 km s^^ kpc~^ and Rqr = 92'.'2 ± 9'.'0 for region A. The 
resonance positions, however, have similar values within the errors. 
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OLR, we hypothesize that the "Q effect" is related to the dynamics of the disk. 

NGC 4939 : [Figures [20] - [21] From the location of the dust lanes in this region, it is 
hard to determine the age gradient direction and hence whether it occurs inside or outside 
corotation. If one assumes it is outside corotation (inverse gradient), the spiral endpoints are 
located near the OLR, as shown in Fig. [201 If, contrariwise, one assumes the age gradient is 
inside corotation, the spiral endpoints may coincide with the location of the 4:1 resonance 
(see Fig. [2Ti) . Region NGC 4939 A has a shape that differs from the models, a fact that 
could be explained if the galaxy metallicity is different from the one assumed (see § 15.31) . 

NGC 3938 : [Figure [22] The signal-to-noise ratio for the optical images of this galaxy 
is lower than for other objects in the sample. The dynamical parameters inferred from the 
data yield resonance positions that do not match the location of the spiral endpoints. 

NGC 7126 : [Figure [23] The inferred locations of major resonances do not match any 
spiral endpoint in the Kg mosaic, which has a low signal-to-noise ratio. However, in the 
i-band mosaic two well defined spiral arms can be seen that extend beyond the position of 
the OLR. We conclude that the gradient featured in region NGC 7126 A may not be due to 
star formation triggered by the spiral shock. 

NGC 1417 : [Figures [21- [26] Regions A, B, and C for NGC 1417 give the same corotation 
position (8 kpc) within the errors (see Table [1]). The spiral arms end at the OLR. 

NGC 7753 : [Figure [27] Apparently, the gradient studied in this object coincides with a 
dust lane feature. Nevertheless, according to the computed the spiral endpoints match 
the OLR within the errors. 

NGC 6951 : [Figure [28] Neither of the resonance positions computed from the data of 
region NGC 6951 A matches the spiral endpoints. 

NGC 5371 : [Figures [29] - [30] The dynamic parameters derived for regions A and B 
result in different locations of the major resonances. Those for region A do not coincide with 
the spiral endpoints of the object. However, region B yields an OLR position that matches 
the spiral endpoints. 

NGC 3162 : [Figures [31] - [32] The computed OLR for region B coincides with the the 
spiral endpoints. The computed errors of the dynamic parameters for region A are larger 
than for region B. 

NGC 1421 : [Figures [33] - [35] Region A has an offset in Q of +0.1 mag when compared 
to regions B and C. On the other hand, regions A and C give resonance positions that agree 
with the spiral endpoints; the best match is for region C, which is also an inverse azimuthal 
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color gradient (i.e., situated after corotation). The feature studied in region B must not be 
caused by star formation linked to the dynamics of the disk. 

NGC 7125 : [Figures[5B]-[55] Region B displays an inverse gradient. The results obtained 
from the three studied regions agree within their errors. The spiral endpoints coincide with 
the OLR in this object. 

NGC 918 : [Figure [39] The stellar population model stretched to the data from region 
A indicates that the spiral arms end at the OLR in this galaxy. 

NGC 578 : [Figures Ho] - HI] Regions A and B, located in two different arms of the four 
that conform this object, yield the same value for the spiral pattern speed. In the deprojected 
mosaics, the northern arms, including the one harboring region A, seem to extend beyond 
the corotation radius without reaching the OLR within the errors. The other arms, where 
region B is located, end at corotation. This is the only object in the sample that presents 
spiral arms ending at corotation. 



Table 4. Observed and derived dynamic parameters. 



Region number Galaxy and region Rn 



(arcsec) 



-Rn 



(kpc) Rend (arcsec) -Rend (kpc) Qoffset (mag) Clp (km s ^ kpc ^) Rck (arcsec) 



RcR (kpc) 



1 

2 

3 

4 

5 

6 
7 
8 

9 

10 

11 
12 

13 
14 

15 
16 
17 

18 
19 
20 

21 

22 
23 



NGC4254 A* 
NGC4254 B 

NGC4939 A* 

NGC3938 A* 

NGC7126 A* 

NGC1417 A 
NGC1417 B* 
NGC1417 C 

NGC7753 A* 

NGC6951 A* 

NGC5371 A 
NGC5371 B* 

NGC3162 A 
NGC3162 B* 

NGC1421 A 
NGC1421 B 
NGC1421 C* 

NGC7125 A 
NGC7125 B* 
NGC7125 C 

NGC918 A* 

NGC578 A* 
NGC578 B 



66.4 ±0.25 
49.3 ±0.25 

129.4 ±0.25 

30.7 ±0.25 
34.6 ±0.15 

22.6 ± 0.25 

25.3 ±0.25 

20.4 ±0.25 

.33.7 ±0.15 

38.8 ± 0.25 

32.2 ±0.25 
41.8 ±0.25 

33.0 ± 0.25 

35.7 ±0.25 

34.6 ±0.15 

57.5 ±0.15 

80.0 ±0.15 

33.6 ± 0.15 

79.1 ± 0.15 

45.3 ±0.15 

32.1 ±0.15 

75.2 ± 0.15 

67.8 ±0.15 



5.3 ±0.4 
3.9 ±0.3 

29.2 ±2.5 

2.3 ±0.2 

7.5 ± 0.6 

6.3 ±0.5 

7.0 ±0.6 
5.7 ±0.5 

11.8 ±1.0 

4.7 ±0.4 

6.8 ±0.6 
8.8 ±0.7 

3.8 ±0.3 

4.1 ±0.4 

4.9 ±0.4 

8.2 ±0.7 
11.4 ± 1.0 

7.3 ±0.6 
17.1 ± 1.5 
9.8 ±0.8 

3.4 ±0.3 

8.3 ±0.7 

7.5 ±0.6 



157.5±7.5 (Ks) 12.6±0.8 



145±5 (Ks) 
100±5 (Ks) 
118.9±2.9 (i) 
60±5 (J) 

95.7±5.8 (i) 
105±5 (Ks) 
115±5 (Ke) 

75±5 (J) 



32.7±2.8 
7.7±0.7 
25.7±2.2 
16.6±1.4 

33.5±2.8 
12.7±1.1 
24.3±2.1 

8.6±0.8 



92.8±2.9 (Ks) 13.2±1.1 

95.7±2.9 (J) 20.7±1.8 

75.4±5.8 {K') 7.9±0.7 

95.7±2.9 (Ks) 10.5±0.9 



-0.13 ±0.05 
-0.16 ±0.05 

0.0 ±0.06 

+0.09 ± 0.06 

0.0 ±0.09 

-0.15 ±0.06 
-0.18 ±0.06 
-0.15 ±0.06 

+0.07 ±0.09 

+0.14 + 0.07 

-0.01 + 0.05 
+0.03 + 0.05 

0.0 + 0.06 
0.0 + 0.06 

+0.05 + 0.09 
-0.04 + 0.09 
-0.05 + 0.09 

0.0 + 0.09 
0.0 + 0.09 
-0.02 + 0.09 

+0.07 + 0.05 

-0.2 + 0.09 
-0.16 + 0.09 



33.5 + 8.2 

47.7 + 11.4 

12.3 + 0.8 

33.8 + 12.9 
19.7 + 2.2 

35.1 + 4.7 

29.9 + 4.1 

33.6 + 5.2 

11.7 + 2.8 
49.9 + 12.8 

31.2 + 7.6 

20.8 + 6.5 

14.7 + 8.3 
25.0 + 8.3 

16.6 + 3.5 
0.5 + 2.2 
22.3+1.7 

9.6 + 4.3 

10.8 ± 1.7 
9.8 + 3.2 

33.4 ±5.7 

15.0 ± 1.9 
15.2 ± 2.2 



to 



81.0 ±8.6 

57.0 ± 5.6 

87.5 ±5.3 
36.9 ±5.2 
40.2 ±4.1 

26.8 ± 2.9 

31.5 ±3.6 

28.1 ±3.6 

51.7 ±9.0 

43.2 ±5.0 

44.3 ±6.1 

66.6 ± 13.3 

81.8 ±33.5 
48.0 ± 7.9 



72.6 ± 13.9 
2230.1 ± 15820.0 

54.0 ±3.6 

69.4 ±22.1 

61.7 ±4.8 

68.5 ± 12.6 

39.6 ±4.9 

92.2 ± 9.7 

90.8 ± 10.9 



6.5 ± 0.7 

4.6 ± 0.4 

19.7 ± 1.2 
2.8 ±0.4 

8.7 ±0.9 

7.4 ±0.8 

8.7 ± 1.0 

7.8 ± 1.0 

18.1 ± 3.2 

5.2 + 0.6 

9.3 + 1.3 
14.0 + 2.8 

9.4 + 3.8 

5.5 + 0.9 

10.3 + 2.0 
316.8 + 2247.2 
7.7 + 0.5 

15.0 + 4.8 
13.3 + 1.0 

14.8 + 2.7 

4.2 + 0.5 

10.1 ± 1.1 
10.0 + 1.2 
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5. Expected Q{rJgi) index profiles in a spiral shock scenario. 

In this section, we consider in a qualitative way a more realistic situation within the 
density wave scenario. 

The relative density distributions of young and old stars after the shock, as well as 
non-circular stellar velocities produce changes in the azimuthal light profiles with respect to 
the idealized case discussed earlier. Also, since we are dealing with different generations of 
stars, metallicity variations should be considered. In what follows, we discuss each one of 
these issues at moderate length, and make an estimate of their impact on the disk dynamic 
parameters derived through the comparison of model and observed Q profiles. 



5.1. Post-shock density and velocity distributions. 



With the purpose of increasing our understanding of gas dynamics in the presence of 
spiral density waves, stati onary spiral shock patterns have bee n studied with both semi- 



analy tical approaches (e.g.. [Roberts 



Slvz et al.ll200a iMartos et al. 



2004 



969l: iGittins fc Clarkell2004r) and numerical simulations 



Yanez et al.ll2008l ). The post-shock density and 



(e.g. 

velocity profiles obtained from such studies show that they depend on many physical param- 
eters. It is reasonable to assume, however, that the newborn stars product of these shocks 
have densities and velocities that are similar to those of the collapsed gas clouds where they 
form, at least in the early stages of their evolution, and even if only a few percent of the gas 
will form stars. Also, it is commonly accepted that dust lanes trace the location of spiral 
shocks. As newborn stars move away from their birth site, different distances are reached 
due to accelerated movements. In figure [5] we show th e gas velocity parallel to the spiral 
equipotential curve from the semi-analytical solution of [Roberta (119691 ): this solution corre- 



sponds to a radius of 10 kpc; a pattern speed, fip, of 12.5 km kpc~ , an arm pitch angle 
i = 8?13; a spiral field strength of 5% that of the axisymmetric field; and a mean gaseous 
dispersion of 10 km s"^. The distance is measured relative to the spiral shock location, and 
the velocity, in the non-inertial frame of reference of the spiral pattern. The variable velocity 
in an inertial frame of reference goes from 235 to 255 km s~^. For the time being, we assume 
that the circular gas velocity is approximately equal to this velocity, and investigate the 
effects of such variable stellar velocity on the 1-D profiles of the photometric index Q{rJgi). 

Figure [6] shows theoretical profiles for the Q{rJgi) index, obtained with the variable 
velocity plotted if Fig. O The models have an IMF upper mass limit of 10 Mq, and a 
constant fraction of young stars of 2% throughout. A comparison with constant velocities of 
120 km s~^ and 240 km s~^ is also shown. 



Table 4 — Continued 



Region number Galaxy and region i?mean (arcsec) Rmean (kpc) -Rend (arcsec) -Rend (kpc) Qoffset (mag) Clp (km s ^ kpc ^) -RcR (arcsec) -RcR (kpc) 



Note. — Col. (3) and (4). Mean radius of the studied regions, in arcsec and kpc, respectively. Col. (5) and (6). Radius of the spiral endpoint, in arcsec and kpc, respectively, and 
bandpass used to determine it. Col. (7). Q index magnitude offset, applied to the solar metallicity model with 2% of young stars in order to fit the data. Col. (8). Spiral pattern speed. 
Col. (9) and (10). Corotation radius, in arcsec and kpc, respectively. 
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Fig. 4. — Ratios of resonance positions, Rres, to spiral endpoints, Rend- Except in the case 
of NGC 4254 (M 99), objects are organized by their Hubble type; values of the resonances 
and spiral endpoints are listed in Table HI Filled circles: OLR; filled squares: corotation; 
empty triangles: 4:1 resonance. 
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Fig. 5. — Gas circular velocity vs. angular distance obtained from the solution by iRoberts 
Jl969h . The distance is measured relative to the spiral shock location, and the velocity, 
relative to the non-inertial frame of reference of the spiral pattern. We assume that the gas 
circular velocity is approximately equal to the velocity parallel to the spiral equipotential. 
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Fig. 6. — Theoretical profiles for the Q{rJgi) photometric index, CB07 models. Solid line: 
stars have a variable orbital velocity, as shown in Fig. [5l dotted line: stars have a constant 
circular velocity of 120 km s~^; dashed line: stars with a constant orbital velocity of 240 km 
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Next, we investigate the effect of including variable, more realistic, mass fractions of old 
and young stars. To this e nd, we use the relative gas density shown in Fig. [3 This density was 



derived by iRobertsI (Il969l ). using the same parameters listed above for his velocity solution. 
We then assume that the fraction of young stars must be 2% at an age of 3 x 10'' years 
(6* ~ 20°), and propagate this fraction to other locations, supposing that the young stars 
share the gas density distribution (for example, the young stars fraction would be about 2.5% 
at an age of 1.5 x 10^ years or 6' ~ 10°). As already stated, for the old stars mass fraction 
we use Pii = 1 - Pj (see equation [H]) . The theoretical Q profiles involving variable stellar 
velocities, and both variable and constant stellar densities are shown in figures [8] and [9l for 
IMF upper mass limits of 10 and 100 Mq, respectively. 

A further possible refinement to the models concerns the rel ative location of the shock 



and the potential minimum. According to lGittins &: Clarkd (120041 ). the shock location moves 



to different azimuthal values for tightly wound spirals. At small radii (inside corotation), the 
shock occurs near the potential minimum; at larger radii, though, the shock weakens and 
moves upstream towards the potential maximum. In real galaxies, the ma ximum surface 



densi ty of old stars and the potential minimum will not coincide exactly (IZhang fc Buta 



20071 ). The gas responds to the potential minimum, while the maximum observed surface 
density of old stars marks the peak of the density wave. The old stars' surface density is 
commonly inferred from obse rvations at 2Atm, al though red supergiants may contribute 20% 



of the flux in this wavelength (IRix fc Riekelll993h . Assuming that the onset of star formation 



occurs almost immediately after the shock, the resulting total stellar density (i.e., considering 
both young and old stars) will be affected by the relative positions of each component. In 
Fig. [TOl we show a possible shape for a density wave taken from the Kg data of NGC 7125, 
assuming all the emission comes from an old population with a constant mass-to-light ratio. 
We try three positions of the density wave peak at, respectively, 10, 15, and 20 degrees away 
from the shock, with increasing widths. With these density distributions of old stars, and 
the variable velocity and relative density distributions for young stars discussed previously 
(the fraction of young stars is taken to be 2% at an age of 3 x 10'' years, and propagated 
to other positions following the gas density distribution shown in Fig. [7]), we produce once 
again theoretical profiles for the Q{rJgi) index. These are shown in figures [TT] and [T^ for 
IMF upper mass limits of 10 and 100 Mq, respectively. The comparison between the more 
complex models and those with constant stellar densities and velocities are also shown in 
this figures, for a fraction of young stars of 2% by mass. From this analysis it is clearly seen 
that galaxy azimuthal light profiles may suffer deformations depending on the density and 
velocity distributions of the underlying stars, being the deformation due to density the most 
important one. 



It is worth noticing that in models with constant stellar density and velocity a degeneracy 
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Fig. 7. — Gas relative density from the shock solution of iRobertd (119691 ). Sq is the unper- 
turbed density in the axisymmetric potential, and Si is the perturbed density due to the 
spiral gravitational field. The resulting total gas density is S = Eq + Si. 
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Fig. 8. — Model profiles for the photometic index Q{rJgi), low IMF upper mass limit. Solid 
line: variable stellar velocity and density as shown in figures [5] and [3, respectively. Dotted 
line: variable stellar velocity and constant density (2% of young stars by mass). The IMF 
lower and upper mass limits are 0.1 and 10 Mq, respectively, for both models. 
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Fig. 9. — Model profiles for the photometic index Q{rJgi), high IMF upper mass limit. 
Solid line: variable stellar velocity and variable stellar density; dotted line: variable stellar 
velocity and constant density. The IMF lower and upper mass hmits are 0.1 and 100 Mq, 
respectively, for both models. 
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0(deg) 

Fig. 10. — Density profiles of old stars for different locations of the density wave relative to 
the spiral shock. Solid line: old population density peak at 10°; dotted line: density peak 
at 15°; dashed line: density peak at 20°. Eq is the unperturbed density in an axisymmetric 
potential; Si corresponds to the density wave perturbation; E = Eq + Ei. Angular locations 
are the same as in previous figures. 
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Fig. 11. — Theoretical Q profiles for variable stellar densities and velocities, in the case of 
an IMF upper mass limit of 10 Mq. Each model corresponds to a different location of the 
stellar density wave peak, as shown in Fig. [TOl Lines have same meaning as in Fig. [TOl except 
for the long- dashed line: constant stellar density and velocity. 
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Fig. 12. — Same as Fig. [11], for an IMF upper mass limit of 100 Mq. The main difference 
witfi respect to Fig. [11] is tlie sliarp peak at ^ ~ 0°. 
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occurs between IMF upper mass limit and fraction of young stars (GG96), where models 
with low IMF upper mass limit and higher young star fraction resemble those with high 
upper mass limit and lower young star fraction (1 %). This degeneracy is shown in Fig. [T3i 
Pi In models with a fraction of young stars that varies with azimuthal position, however, 
this degeneracy is broken, as is shown in figures [TT] and [T2l The models with an upper mass 
limit of 100 Mq show a peak close to the shock that is absent in those with Mapper = 10 Mq, 
although it is true that the peak might be lost in data with a poor signal-to-noise ratio. 



5.2. Non-circular motions. 



The motions of young stars under t he ass u mption o f spiral density wa v e triggering have 



been st udied by, amo ng others, lYuanl (119691 ): IWielenI (Il973l . Il978l . Il979l ): [Fernandez et al. 



( I2OO8I ). IWieleru (119791 ) has emphasized that, in a galaxy with spiral density waves, newly-born 
stars have not had time to reach dynamical equilibrium with the galaxy potential. Hence, 
instead of responding to a stationary wave like the old population, young stars migrate out 
of the arms following complicated, non-circular, orbits. The aforementioned trajectories are 
not necessarily closed; also, they can run almost parallel to the arms for significant stretches, 
with the consequence that stars are not seen leaving the arms until they have somewhat aged 
and at a radius where the dust lane location does not mark the site of star formation; stars 
might even move to locations upstream of the shock. The color gradients discussed before 
would overlap with stars drifting back to the arms. 

Another consequence in real galaxies is that orbits are no longer vertical lines in the 6 
vs. In R plane; rather, they are oblique and wavy, as shown in Fig. (TH Given that we have 
collapsed our data in radius to improve the signal-to-noise ratio, our analysis is insensitive 
to this effect. 

Finally, if molecular clouds suffer the effects of the spiral shock, drift velocities will be 
smaller and s tars will take longer to migrate away from the arms, as shown in Fig. 5 of 
WielenI ( 119791 ): there, lines of equal stellar age are closer together than in his Fig. 4, where 
clouds and young stars keep their pre-shock velocities. 

Qualitatively, these phenomena are akin to the variable stellar velocities and densities 
considered in § 15. H and whose impact on the derived density wave parameters is minimal, 
as briefly discussed in § 15.41 Numerical simulations beyond the scope of this paper would be 



A similar degeneracy exists between IMF upper mass limit and length of star formation burst (GG96), 
in the absence of independent constraints on their values. 
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Fig. 13. — Degeneracy in constant stellar density models between IMF upper mass limit 
and fraction of young stars. Solid line: Mupper = lOOM© and 1% (by mass) of young stars; 
dotted line: Mupper — IOMq and 2% of young stars. Burst duration is 2 x 10^ years. The 
difference in the derived from the two models would be ~ 3 km s"^ kpc~^. 
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required to assess their repercussions in a quantitative way. 



5.3. The role of metalhcity. 

The chemical composition of spiral dis ks has been studied by several authors, mainly 

through spectroscop i c studies of HII regio ns Jvila-Costas k Edmunds lll992l:l IZaritskv. Kennicutt fc Huchral 
19941 : Ivan Zee et alJEoOsI : [Camettl EqJ • IZaritskv. Kennicutt fc Huchral Jl994h measured 
the oxygen abundancej^^l using a photoionization model for calibration, in 39 disk galaxies 
at the characteristic radius r = 0.4po, where po is the B = 25 mag arsec"^ isophotal radius. 
Their main result for various Hubble types is shown in Fig. [T3 

Adopting the usual normalization X + Y + Z = 1, where X, Y, and Z are, respectively, 
the abundances per unit mass of hydrogen, helium, and the remaining elements, we have: 



Z{0) 
X 



IQiO/H), 



(10) 



where Z{0)/Z is the fraction of Z due to oxygen. This fraction va ries with metallicity 
between 41% and 53% in the Local Group of galaxies (jPeimbertI l2003l ) . Taking the values 
X = 0.75 and Z{0)/Z = 0.45, the expression for Z in terms of the oxygen abundance is: 



log(Z) ~ 1.43 + \og{0/H). (11) 

Assuming that young stars have metallicity values around the ones inferred from abun- 
dance ratios, we obtain that Z is in the range [0.006, 0.05] for the disk stars in the Zaritsky et 
al. (1994) sample. In Fig. [TBlwe show the behavior of the Q index for various combinations 
of metallicities for the young and old stars, using the CB07 models. In these models the 
adopted solar metallicity is Zq = 0.02. The duration of the burst is 2 x 10^ years, both 
young and old stars have a Salpeter IMF with Miower = O.IM© and Mupper = lOOM©. The 
fraction of young stars is 2% by mass, and the old stellar background population is 5 x 10^ 
years old. 

When young stars have Z > Zq, Q reaches much higher values than in models where 
young stars have subsolar metallicities. Models that differ in the metallicities of the young 
population show different Q values mainly between 1x10^ and 3x10'' yr, when the young 



-"^^Defined as 12 + log(0/H), where (0/H) is the number ratio of oxygen to hydro gen atoms. The present 
day local ISM has the value 12 + log(0/H) = 8.79 ± 0.08 (jCarigi fc Peimbertlljooi) . 
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stars are most prominent; if, on the other hand, the young stars have the same Z but the 
old populations do not, the models run basically parallel, with a roughly constant offset in 
the Q values at all ages. 



5.4. Simple vs. complex models. 

Through the comparison between examples of the more complex models presented and 
our data, we estimate that ignoring the deformations produced by variable densities, variable 
velocities, and different metallicities (see § 15.11 § 15.21 and § 15.31) would translate into a 
maximum error of approximately ±1 km s~^ kpc~^ in fip (see also the Appendix). This 
quantity is well within the random errors computed for f2p shown in Table HI 

The small size of this error might signal a sort of selection bias: if non-linear effects 
were very important, it would be unlikely to find a gradient. On the other hand, remarkably, 
we have been able to detect gradients in an unprecedented 10 out of 13 galaxies. Not only 
have we quadrupled the historic number of detections, but we have found gradients in more 
than 75% of the subsample analyzed in this paper, that does not even comprise (with the 
exception of M 99) the "Q effect" galaxies, i.e., the objects with the strongest features. This 
means that non-linear effects that would deform all gradients beyond recognition might not 
be too significant in galaxies with spiral density waves. When GG96 found the gradient 
in M 99, an obvious question that arose was whether NGC 4254 was an exception or if an 
adequate technique to search for the gradients had been finally devised. The present work 
answers that question unambiguously. 

Furthermore, the comparison of the gradients with the simplest of models yields orbital 
resonance positions that match the observed spiral endpoints. These results agree with the 
predictions of density wave theory, and establish a strong link between disk dynamics and 
large scale star formation, as we argue below. 



Discussion and conclusions. 



If star formation is related to disk dynamics then, according to theory, the Rres/Rend 
ratio must be close to o n e, for either the 4:1 resonance, the corotation radius, or the OLR. 
Contopoulos &: Grosb0ll (Il986l ). based on orbital calculatio ns, determined that the spiral 
patter n of strong spirals must extend to the 4:1 resonance. iPatsis. Contopoulos fc Grosb0l 
(Il99ll ) defined strong spirals as those with large pitch angles (i.e., Sb and Sc galaxies, that are 
not tightly wound) ; their theory is considered "non linear" . Conversely, the "linear" theory 
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of spiral density waves concerns itself with tightly wound sp irals and has concluded that 
the arms of some normal (nonbarred) sp irals reach corotation (|Linlll970f) and are sta tionary. 



tne arms oi some normal (^nonbarrea ) sp irals reacn corotation (|ljinlll^i^ui) and are sta tionary, 
while those of others gro w to the OLR jMarklll97d : llin fc LaulllQTol : lToomrelll98ll ). In the 



Contopoulos fc Grosb0ll Jl986h treatment, the "linear" theory is recovered when spiral 
are not strong (as in Sa galaxies). 



arms 



In Fig. m we can see that most spirals in our sample extend to the OLR and one of 
them (NGC 578) reaches corotation. The mean of the ratio Rqlr/ Rend, for all 13 points, 
is 0.95±0.03. The reduced x^/^ of this result is 7.120 the probability of this result being 
due to chance for 13 degrees of freedom is less than 1 out of 10,000. If we remove the points 
corresponding to NGC 3938, NGC 7126, NGC 6951, and NGC 578, Rqlr/ Rend = 0.98±0.04, 
with a reduced x^/""- = 0.48, for n = 9. This last result may indicate that some of our errors 
are overestimated. It is interesting to note that none of our objects is an Sa (weak spiral) 
galaxy. From this we may conclude that the "linear" result for the extent of spiral patterns 
applies to strong spirals too! 

On the other hand, for most objects the spiral begins at the location of the ILR, as 
expected. For the analysis presented here, though, the location of this resonance should be 
taken with care, since we are employing fiat rotation curves, even at small radii. When using 
real rotation curves, the positions of the ILRs for our sample may vary. The possibility also 
exists that there is no ILR in some objects. 

In some regions, the downstream values of the Q index are lower than the models. 
This is the case for NGC 4939 A, NGC 3162 B, NGC 1421 A, NGC 1421 C, NGC 7125 A, 
NGC 7125 C, NGC 578 A, and NGC 578 B. For the present analysis we have assumed that 
every star-forming region maintains its previous circular motion after the spiral shock. As 
already stated, in real galaxies the situation is different: variations in the v elocity vector are 
present due to the shock itse lf and the resulting loss of angular momentum (lYuan fc Grosb0l 
198ll : [Fernandez et al.ll2008l ). This effect could explain the "downstream fall" of the gradi- 
ents. Metallicity is another factor that can play a significant part: a lower metallicity of the 
older stars may produce a steep drop of the gradients (see § 15.31) . 

We have shown that azimuthal color gradients are common in spiral arms of disk galax- 
ies. Even inverse color gradients have been found in the arms of NGC 1421, NGC 4939, 
and NGC 7125. The observed picture is, of course, not as clean as originally envisioned, 
for several reasons. Dust and HII regions may mask the gradients (GG96); indeed, all the 
detected gradients have been satisfactorily fitted with models where Mapper = 10 Mq. Other 
star formation mechanisms, like self-propagating star formation, may take place simultane- 



For an expected value of Rqlr/ Rend — 1- 
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ously with density wave triggering in spiral arms and disks in general. Likewise, there may 
be substructure formatio n, as a result of non-liii e ar and chaotic effects associated to the 
wave phenomenon itself (jChakrabarti et al.l l2003l : iDobbs fc Bonnelll l2006l : iKim &: Ostriker 
20061 : IShetty &: Ostrikerll2006l ). In the future, three aspects can be improved when carrying 
out this type of study. (1) A better determination of the inclination angle of the galaxy, 
and of the rotation curve at different radii; this will reduce the uncertainty in Qp. (2) The 
inclusion in the models of the variable density distributions of young and old stars, as well 
as of the stellar velocity changes near the spiral shocks. This will account for non-circular 
motions near the spiral shocks and for the higher order terms in equation [2l Numerical 
simulations or semi-analytical treatments may be needed, particularly for the young stars; 
also, the contribution of red supergiants to the observed i^'-band surface brightness should 
be considered when modeling the distribution of old stars. (3) Spectroscopic studies could 
be very helpful as diagnostics of the population properties, particularly of the metallicity 
along the gradients. 

According to the dynamic parameters derived from our analysis, spiral arms mostly 
extend to the OLR and sometimes to the corotation radius. Ten of thirteen objects (77%) 
support this conclusion. Similar results ha ve been obtained by other authors as well (e.g., 
Elmegreen et al.lll992l : IZhang fc Butall2007l ). In view of the consistency between the pattern 
speeds and resonance positions determined from the azimuthal stellar gradients, and the 
predictions of density wave theory, we conclude that disk dynamics do play an important 
role in large scale star formation in some spiral galaxies. 
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A. Pattern speed and corotation radius error calculation. 

In principle, the three contributors to the random error in the pattern speed, a^p (km 
s~^ kpc^^), are the uncertainties in: the inclination angle, cTq,; the rotation velocity, cr^„^^; 
and the distance to the galaxy, ctd. 

In order to compute an^ from equation [6l and hence the uncertainty in the resonance 
positions, we replace Rmean, Vrot, and d with their corresponding expressions in terms of 
independent variables. We obtain: 
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a,,^^-^eos-( ^t"'M . (Al) 
rVisma) t v ' 



(A2) 



where constant ki = ^^^220 (pixels), Ps is the adopted plate scale in arsec pixel"^, constant 
^2 = 9.766 X 10^ (yr km s^^ kpc^^), r is the mean deprojected radius of the region in pixels; 
Vmax is the maximum rotation velocity in km s~^; {xi,yi) and {x2,y2) are the cartesian 
coordinates (in pixels) of the extremes of the curved line segment corresponding to the 
studied region in the non-deprojected image (the origin is located at the center of the object 
and the galaxy major axis runs along the y coordinate); D is the distance to the object in 
kpc; a is the inclination angle used to deproject the image; and t is the stellar model age in 
years. 

The corotation radius, Rqr (kpc), is: 



RcK = j^^. (A3) 
(sm a)\lp 

The expressions for each contributor to the error calculation are: 



AQp[a] = Qp{a + aa, Vmax, D,t) -Q D,t) (A4) 

A^p[Vmax] = ^p(a, Vmax + CTy^^^.D, t) - Qp{a, Vmax, D, t) (A5) 

/^VLp[D] = Qp{a, Vmax, D + an, t) - Qp{a, Vmax, D, t) (A6) 
an, = y^(A^]p[a])2 + {A^plvmaxW + (A^]p[Z^])2. (A7) 

Values for ±a were calculated, and the absolute AQp values were averaged. Similar expres- 
sions hold for (JR(^^. Figure fTTl shows the AQp and Ai^cR absolute average values for each 
region listed in Table |H The systematic error due to the choice of stellar population model 
was estimated by comparing the Qp values derived with the simple models of § [3], on the one 
hand, with examples of the more complex models discussed in § 15. 4| on the other. 
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Fig. 14. — Left: Orbit in cartesian coordinates (arbitrary units). Right: The same orbit in the 
6 vs. InR plane. Vertical lines represent circles every ~ 0.3 arbitray units (R — y^x"^ + y^). 
Horizontal lines are distributed every 90°. 
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Fig. 15. — Characteristic abundanc e vs. Hubble type (T-type as defined by RC2) from 
Zaritsky. Kennicutt fc Hucfiral (Il994j ). Smaller T-types correspond to earlier type galaxies. 
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Fig. 16. — Q index behavior with different combinations of metaUicities for the young and 

0.004; dotted line: Zy^nng = 0.008, Z^ia = 0.004; 



old stars. Solid line: ^young = 0.004, Zoid 



dashed line: Z, 
dashed line: Z, 



young 



young 
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0.004; dashed-dotted line: Zyoung = 0.02, = 0.02; long- 
■ 0.004; long-dashed-dotted line: -^young = 0.05, Z^i^ — 0.02. 
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Fig. 17. — Absolute AQp and Ai?cR average values for galaxy regions. Regions are identified 
by their numbers in Table HI Solid line: contribution from inclination angle, a; dotted line: 
contribution from rotation velocity, Vmax] dashed line: contribution from distance to the 
galaxy, D. The measurements are shown as continuous lines for easiness of reading; in 
reality, they are discrete points. 
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Fig. 18. — Region NGC 4254 A. Panel a: Deprojected spiral galaxy NGC 4254 in the optical 
(7-band. The display is in logarithmic scale. The studied regions are marked on the figure. 
Panel h: 1-D data vs. azimuthal distance for region NGC 4254 A. Solid line: Q index; dotted 
line: [g — J) color; dashed line: i^^-band profile. Panel c: Zoom-in of region, and 1-D Q index 
vs. azimuthal distance. Solid line: data; dotted line: model. Data error bars and direction 
of rotation are indicated. Panel d: i^"s-band deprojected mosaic of spiral galaxy NGC 4254. 
Solid line circles: location of the ILR and the OLR, as obtained from the comparison between 
data and stellar population SPS model shown in panel c. Dashed line circles: ±lcr errors. 
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Fig. 19.— Region NGC 4254 B. 
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Fig. 20. — Region NGC 4939 A. In order to fit the model to the data we assume the gradient 
is inverse, occurring beyond corotation. This can also be seen from the fact that time 
evolution (top axis in panel c) and rotation have opposite directions. 
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Fig. 21. — Region NGC 4939 A. In order to fit the model to the data we assume now that 
the gradient occurs inside corotation. In this case, time evolution (top axis in panel c) and 
rotation have the same direction. Panel ci, solid line circle: location of the 4:1 resonance, as 
obtained from the comparison between data and SPS model shown in panel c. 



- 55 - 




Fig. 22. — Region NGC 3938 A. Panel a: Deprojected spiral galaxy NGC 3938 in the optical 
r-band. 
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Fig. 23. — Region NGC 7126 A. Panel d: i-band deprojected mosaic of spiral galaxy 
NGC 7126. 
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Fig. 24. — Region NGC 1417 A. Panel d: J-band deprojected mosaic of spiral galaxy 
NGC 1417. 
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Fig. 25. — Region NGC 1417 B. Panel d: J-band deprojected mosaic of spiral galaxy 
NGC 1417. 
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Fig. 26. — Region NGC 1417 C. Panel d: J-band deprojected mosaic of spiral galaxy 
NGC 1417. 
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Fig. 27. — Region NGC 7753 A. Panel d: i-band deprojected mosaic of spiral galaxy 
NGC 7753. 
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Fig. 28.— Region NGC 6951 A. 
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Fig. 29.— Region NGC 5371 A. 
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Fig. 30.— Region NGC 5371 B. 
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Fig. 31. — Region NGC 3162 A. Panel d: J-band deprojected mosaic of spiral galaxy 
NGC 3162. 
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Fig. 32. — Region NGC 3162 B. Panel d: J-band deprojected mosaic of spiral galaxy 
NGC 3162. 
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Fig. 33.— Region NGC 1421 A. 
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Fig. 34. — Region NGC 1421 B. Panel d, solid line circle: location of the ILR and the OLR, 
as obtained from the comparison between data and SPS model shown in panel c. The errors 
for the resonance positions are much greater than the computed resonance radii themselves 
(see Fig. [T71) . This is due to the fact that region B lies right on the corotation position 
obtained from regions A & C. 
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Fig. 35.— Region NGC 1421 C. 
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Fig. 36. — Region NGC 7125 A. Panel d: J-band deprojected mosaic of spiral galaxy 
NGC 7125. 
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Fig. 37. — Region NGC 7125 B. Panel d: J-band deprojected mosaic of spiral galaxy 
NGC 7125. 
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Fig. 38. — Region NGC 7125 C. Panel d: J-band deprojected mosaic of spiral galaxy 
NGC 7125. 
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Fig. 39. — Region NGC 918 A. Panel d: K'-band deprojected mosaic of spiral galaxy 
NGC 918. 
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Fig. 40. — Region NGC 578 A. Panel d, solid line circles: location of the ILR and corotation 
radius, as obtained from the comparison between data and SPS model shown in panel c. 
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Fig. 41. — Region NGC 578 B. Panel d, i^'s-band deprojected mosaic of spiral galaxy 
NGC 578. solid line circles: location of the ILR and corotation radius, as obtained from the 
comparison between data and SPS model shown in panel c. 



